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Abstract 

The localization length has been derived for one-dimensional bi-layered structures with random perturbations in the 
refractive indices for each type of layers. Main attention is paid to the comparison between conventional materials 
and those consisting of mixed right-hand and left-hand materials. It is shown that the localization length is described 
by the universal expression for both cases. The analytical results are confirmed by numerical simulations. 
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Introduction - In many fields of fundamental and 
applied physics much attention is paid to the wave prop- 
agation and electron transport in one-dimensional peri- 
odic structures with elementary cells consisting of two 
components (see, e.g. OJJ] and references therein). For 
example, it could be two optical or electromagnetic ma- 
terials, or the pair of quantum wells and barriers in elec- 
tronics. The interest to such bi-layer structures is due to 
various applications for which one needs to know how 
to create materials, metamaterials, or semiconductor su- 
perlattices with given transmission properties. 

One of the important problems that still remains 
open, is the role of a disorder that cannot be avoided 
in experimental devices due to fluctuations of the width 
of layers or due to variations of the medium parame- 
ters, such as the dielectric constant, the magnetic per- 
meability, or the barrier hight for electrons. In spite of 
remarkable progress in this field, the majority of stud- 
ies of the wave (electron) propagation through random 
structures are based on various numerical methods, with 
an additional assumption of rapidly deca ying correla- 
tions Q 1 fl S SB S H M M El HE As for 
the analytical results, they are mainly obtained either for 
systems with a white-noise disorder 11411 . or for the pat- 
terns with correlated disorder, however, with delta-like 
potential wells [15] or barriers |16U17ll . 



In this paper we derive the unique analytical expres- 
sion for the localization length L; oe that is valid for a 
discrete bi-layer structures with weakly disordered re- 
fractive indices of both basic slabs. In contrast with the 
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general approach, see 111 811 . here we are interested in a 
particular problem of comparison of conventional pho- 
tonic crystals with metamaterials. 

Model - We consider the propagation of an electro- 
magnetic wave of frequency u> through an infinite array 
(stack) of two alternating a and b layers (slabs). Every 
kind of slabs is respectively specified by the dielectric 
permittivity s a magnetic permeability jj. a j,, refractive 
index n a , b = ^e a ^ a ,b, impedance Z fl , ft = n a ,bln a ,b and 
wave number k a ^ = a>n a j,lc. We address two cases 
when a slabs contain conventional right-handed (RH) 
optic material, while b layers are composed of either 
RH or left-handed (LH) material. The combination of 
RH-RH slabs is caled homogeneous stack whereas the 
array of RH-LH layers is called mixed stack. Following 
Ref. Illll . a disorder is incorporated via the dielectric 
constants only (compositional disorder), so that 

[X a = 1, n a (n) = 1 + rjjn); (la) 
fX h = ±1, n b (n) = ±[1 + rjb(n)]. (lb) 

Here integer n enumerates the elementary (ab) cells. 
The upper sign is related to RH material while the lower 
one is associated with LH media. Every alternating slab 
has the constant width d a or db, respectively. 

Without disorder, r] a ,b{n) = 0, all layers are perfectly 
matched (Z a = Zb) and the stack is equivalent to the 
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homogeneous medium with the refractive index n, 



con/c, 



n-\d a ± d b \/(d a +d b ), 



(2) 



with no gaps in the linear spectrum. Remarkably, in 
the ideal mixed stack (e fl = \i a = 1, s b — fi b = — 1, 
Z a = Z b = 1) with equal slab widths, the phase velocity 
c/n diverges. Therefore, our consideration is meaning- 
ful only when two layers have different widths d a + d b . 

The random sequences rj a (n) and rj b (n) are statisti- 
cally homogeneous with the zero mean, (j] a j,{n)) = 0, 
and binary correlation functions defined by 

(r] a (n)T] a (n')) = a 2 , K a (n - «') , (3a) 
(r]b(n)T] b (n')) = oj, K,,(n - n') , (3b) 
(Va(n)rib(ri')) = a 2 ^ K ab (n - n). (3c) 

The average (...) is performed over the whole array or 
due to the ensemble averaging, that is assumed to be 
the same. The auto-correlators K a {r) and K b (r) as well 
as the inter-correlator K ab {r) are normalized to one: 
K a (0) = K b (0) = K ab (0) = 1. The variances o- 2 a and 
cr 2 b are obviously positive, however, the term <j 2 ab can be 
of arbitrary value. We assume the disorder to be weak, 
(k ab d a 

,b°"a,b) 2 1, that allows us to develop a proper 
perturbation theory. In this case all transport properties 
are entirely determined by the randomness power spec- 
tra 'Ka(K), 'Kb(k), and < K ab (k), defined by the relations 

DO DO 

<K(k) = ^ K(r) exp(-ikr) = 1 + 2 V K(r) cos(fcr). 

r=-oo r=l 

(4) 

By definition (0, all the correlators K a (r), K b (r) and 
K ab (r) are real and even functions of the difference 
r = n — n' between cell indices. Because of this fact 
and due to their positive normalization, the correspond- 
ing Fourier transforms < K a (k), "Khik) and < K ab (k) are real, 
even and non-negative functions of the dimensionless 
lengthwise wave number k. 

Within every a or b layer, the electric field of the wave 
obeys the ID Helmholtz equation with two boundary 
conditions at the interfaces between neighboring slabs, 



-[£+<,bWa,b{x) = Q, (5a) 

l/taiXi) = fib(Xi), l*a VlOi) = Mb Vfcfe)- (5b) 

The x axis is directed along the array of bi-layers, and 
x = Xi stands for the interface coordinate. 

Hamiltonian map approach |l5l EH H - The so- 
lution of Eq. © can be presented as the recurrent rela- 
tions for the wave function \jf m = Q„ and its derivative 



(c / a>)t//' an = P„ at the two opposite edges of the nth ele- 
mentary (a, b) cell, 

Gn+i = A n Q„ + B„P„, P„+\ = -C n Q„ + D n P n . (6) 

The factors A„, B„, C„, D n read 

A„ = cos Tp a cos Tp b - Z~ x Z b sin^ a sin <^,, (7a) 

B„ = Z a sin Tp a cos Tp b + Z b cos Tp a sin Tp b , (7b) 

C„ = Z^ 1 sin^ fl cos^/, + Z h 1 cos^ fl sin^,, (7c) 

D„ = cos Tp a cos Tp b - Z a Z h y sin lp a sin^. (7d) 

They depend on the cell index n due to the random re- 
fractive indices ([TJ, which enter the impedances Z a tb , as 
well as due to the phase shifts 



<Pa,b(n) = f>a),W + T} a ,b(n)], 

tfia = i0d a /c, tp b = +LL>d b /c. 



(8) 



It is noteworthy to emphasize that the recurrent relations 
(O can be treated as the Hamiltonian map of trajectories 
in the phase space (Q, P) with discrete time n for a linear 
oscillator subjected to time-depended parametric force. 

Without disorder r\ a ,b(n) = 0, the factors ^ do not 
depend on the time n. Therefore, the trajectory Q„,P n 
creates a circle in the phase space (Q, P) that is an image 
of the unperturbed motion, 



Q,,+i = Qn cos y + P„ sin y, 
P n +\ = -Q n sin y + P„ cosy. 



(9) 



The unperturbed phase shift y over a single elementary 
(ab) cell is defined as 



y = <p a + <p b = a>(d a + d b )/c, 



(10) 



that gives y = K(d a + d b ) due to Eq. Having the 
circle, it is suitable to pass to action-angle variables R n 
and 6„ via the standard transformation 



Q n = R„cos6„, 



P n — R„ sin 6„ . 



(11) 



By direct substitution of Eq. (fTTT > into the map (0, one 
can reveal that for the unperturbed trajectory the radius 
R n is conserved, while its phase 6 n changes by the Bloch 
phase y in one step of time n, 



Rn+i - Rn, 



"n+l — v n -y. 



(12) 



Evidently, a weak random perturbation results in a 
small distortion of the circle ( fT2b that can be evaluated 
in the following way. First, in the initial map (O we ex- 
pand the factors (0 up to the second order in the pertur- 
bation parameters rj a , b (n) <K 1 entering the impedances 



Z a b(n) and phase shifts Tp a Jji). After getting the per- 
turbed map for Q„ and P„, we pass to action-angle vari- 
ables R„ and 6„ with the use of Eq. (fTTT i. All these 
quite cumbersome calculations allow us to derive the 
perturbed map for the radius R„ and angle 9„ keeping 
linear and quadratic terms in the perturbation: 

R l + il R l = 1 + 1a(n)V a (0„) + m(n)V b (e n ) 
+jfynWa + nl(n)W b + 7] a (n)71b(n)W ab , (13a) 
6„+i -0 n +7 = Va{n)U a {e n ) + rib(ri)Ub(e n ). (13b) 

Here the functions standing with random variables 
Tja,b(n) are described by the expressions: 

V a (9 n ) = -2 sin (f a sin(2#„ - ip a ), (14a) 

V b (9„) = -2 sin (pi, sm(26„ - y - <p a ), (14b) 

W a = 2 sin 2 <p a , W h = 2 sin 2 <p b , ( 1 4c) 

W„b = 4 sin (p a sin <p b cos y; (14d) 

U a (6„) = - sin <p a cos(2fl„ - y> a ), (14e) 

Ub(0„) = - sin 5p 4 cos(26>„ -y-<p a ). (14f) 

Note that in Eqs. ( TT4T > we keep only the terms that con- 
tribute to the localization length L/ oc in the first non- 
vanishing order of approximation. Since in Eq. ( 1 1 3ab 
the factors V„j, containing 8„ are always multiplied by 
T] a ,b(n), only linear terms in these perturbation param- 
eters are needed in the complementing recurrent rela- 
tion (1 1 3bb for the angle 0„ . Relations ( [TBI constitute the 
complete set of equations in order to derive the localiza- 
tion length of the system under consideration. 

Localization length - We define the localization 
length L[ oc via the Lyapunov exponent A |T^, 20ll . 



d a + d h 



-'he 



The average {ab) is performed over the disorder param- 
eters T] a ^{n), and the average ab is carried out over the 
random phase 9„. Now we substitute the recurrent rela- 
tion ( 1 1 3ab into the definition (TT~5T > and expand the loga- 
rithm within the quadratic approximation in the pertur- 
bation parameters. Then we perform the averaging over 
both the disorder and rapid phase assuming the distri- 
bution of 6 n to be homogenous within the first order of 
approximation. After some algebra we arrive at the final 
expression for the Lyapunov exponent: 



A = -a 2 c p< a {2y) sin 2 <p a + -cr 2 b %,(2y) sin 2 <p b 



+ vib K ab(2y) sin ip a sin <p b cosy. (16) 

Note that the Eq. ( TToT ) is expectedly symmetric with re- 
spect to the permutation of slab indices a <-> b. 



In accordance with Eq. ( [ToT l the Lyapunov expo- 
nent A (the inverse localization length L,~J,) consists 
of three terms. The first two terms are contributed 
respectively by the correlations between solely a or 
solely b slabs. Therefore, these terms contain the auto- 
correlators < K a (2y) or r K b (2y). The third term emerges 
due to the inter-correlations between two, a and b, dis- 
orders. It includes the inter-correlator < K a h(2y). 

Eq. (fT5T > manifests that the only difference for the ho- 
mogeneous RH-RH and mixed RH-LH stacks is due to 
the sign in the phase shift <p b . This affects the value 
(TTOl i of the Bloch phase y and the sign at the third inter- 
correlation term. One can also see that the Lyapunov ex- 
ponent typically obeys the conventional frequency de- 
pendence 



A oc L h l c cc co 2 



when 



0. 



(17) 



However, specific correlations in the potential, taken 
into account in Eq. ( TTSI l. may result in a quite unusual 
^-dependence, see B18I1 . 

The Lyapunov exponent A(a>) exhibits the Fabry- 
Perot resonances associated with multiple reflections in- 
side a or b slabs from the interfaces. As is known, they 
appear when the width d a , b of corresponding a or b layer 
equals to an integer multiple of half of the wavelength 
2ncju) inside the layer, 

cd/c - s a njd a or a>/c = s b n/d b , s aib = 1,2,3, 

(18) 

At the resonances the factor sin ip a or sin tpb in Eq. ( [ToT l 
vanishes, resulting in the resonance increase of the lo- 
calization length L[ oc and consequently, in suppression 
of the localization. In the special case when the ratio of 
slab widths d a and dh turns out to be a rational number, 
d a jdb = Sal Sb, some resonances from different types of 
layers coincide and give rise to the divergence of the lo- 
calization length Li oc {to). Remarkably, the Fabry-Perot 
resonance is quite broad because it is caused by vanish- 
ing of smooth trigonometric functions i2lll . 

Of special interest are long-range correlations leading 
to the divergence or significant decrease of the localiza- 
tion length Li„ c (u>) in the controlled frequency window. 
This effect is similar to that found in more simple ID 
models with correlated disorder. In our model this effect 
is due to a possibility to have the vanishing values of all 
Fourier transforms, = Kb — "Kab = 0, in some inter- 
vals of frequency a>. For example, one can artificially 
construct an array of random bi-layers with such power 
spectra that abruptly vanish within a prescribed inter- 
val of co, resulting in the divergence of the localization 
length (l5 , 22 , 3 12 ] • On the contrary, with the use of 



specific correlations one can decrease the localization 
length, and significantly enhance the localization 023I1 . 

Finally, let us consider the particular case of the 
white-noise disorders for a and b slabs, 

c^ = 0, <K a (k) = <K b (k) = L (19) 

Here, the Lyapunov exponent and the inverse localiza- 
tion length turn out to be exactly the same for both ho- 
mogeneous RH-RH and mixed RH-LH stack-structures, 

d a + db 1/2-2 2-2 \ ,, m 

A = — = -(erf, sin 2 <p a + (r l b sin 2 <p b ). (20) 

The numerical results shown in Fig. 1 are obtained with 
the use of Eqs. ©, without any approximation. In the 





Figure 1 : Lyapunov exponent versus frequency. Top: RH-RH (circles) 
and RH-LH (triangles) media for <r a ss <Tt ~ 0.006 , d a = 0.6 , d\, = 
0.4,c = 1 and the length of sequence is N = 10 6 . Bottom: the same 
for cr a x crb ~ 0.3 and the N = 100. Smooth curve depicts Eq. 120) . 



higher panel one can see that for a very long sample 
and weak disorder the analytical expression ( f20b per- 
fectly corresponds to the data, apart from fluctuations. 
For each case only one realization of the disorder was 
used. The fluctuations can be smoothed out by an addi- 
tional ensemble averaging. In order to see whether our 
predictions can be used in experiment, we also show in 



the lower panel the data for a very short sample and very 
strong disorder. As one can see, the analytical result is 
also valid for small frequencies, and gives the qualita- 
tively correct Lyapunov exponent for large values of to. 

The more detailed comparison for a>(d a + db)/c <K 1 
and a>(d a + db)/c » 1 also shows a nice correspon- 
dence. In this respect, a very unusual result obtained in 
Ref. fllll . namely, A ~ co 6 , seems to be entirely related 
to a specific case of equal widths, d a — db- 
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